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We show that space-time evolution of one-dimensional fermionic systems is described by nonlinear 
equations of soliton theory. We identify a space-time dependence of a matrix element of fermionic 
systems related to the Orthogonality Catastrophe or boundary states with the r-function of the 
modified KP-hierarchy. The established relation allows to apply the apparatus of soliton theory to 
the study of non-linear aspects of quantum dynamics. We also describe a bosonization in momentum 
\Q ■ space - a representation of a fermion operator by a Bose field in the presence of a boundary state. 

o 
o 

(N ■ I. INTRODUCTION 



o 



In the seminal paper by Date, Jimbo, Kashiwara and Miwa Q] Sato's approach to integrable hierarchies of soliton 
equations has been formulated in terms of auxiliary fermions. In this formulation the hidden symmetries of soliton 
equations become explicit symmetries of free fermions, the modes of which are labelled by one-dimensional parameter 
with respect to the algebra of fermionic bilinear forms - gl(oo). It was shown that the r-function of soliton theory can 
be expressed as a certain matrix element of lD-fermions. The paper |]J and subsequent works of the Kyoto school 
(see @ and references therein) laid out a mathematical foundation of bosonization — a commonly used representation 
of fermionic modes in terms of a Bose field. These works established a deep relation between quantum physics and 
c"| ' classical nonlinear equations. 

In this paper we adopt this theory to a particularly interesting area of electronic physics - dynamics of fermionic 
coherent states propagating in real space and time generated by the Galilean Hamiltonian 
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We will show that the real time evolution of important objects of electronic physics also obeys classical non-linear 
soliton equations. Therefore the powerful apparatus of soliton theory can be used to study complex aspects of non- 
equilibrium physics of effectively one-dimensional quantum systems. We will report some applications of the results 
of this paper to non-linear electronic transport elsewhere. Here we focus on the formal connection of soliton theory 
and quantum dynamics. In this course we encounter an alternative bosonization with respect to the boundary state, 
but now in momentum space. 

Before formulating our main results, we introduce notations and give a brief review of some important facts on 
a connection between solutions of soliton equations and fermionic matrix elements. The literature on the subject is 
enormous, we restrict the reference list to the original papers 0,0- 

X 

» ! ■ A. Bilinear identity 

a | 

Consider a Fermi-system on a unit circle or, equivalently, ID fermions subject to periodic boundary conditions. 
We let ip p be a fermionic mode of the fermionic field ip(x) — J2 P e^ipp, labelled by integer momenta, p. The 
conjugated field is given by ip^x) = Ylp e ~ ipx,, fip with basic (anti-) commutation relations {tp P ,tp q } = {vJp,Vj\} = 
and {ip p , ?/>J} = S PtQ . The ground state of the system is a filled Fermi-sea such that all modes in the interval between 
the Fermi-points are filled. We will be interested only in one, say, the right chiral sector where all essential modes 
are concentrated around the Fermi point +pf- It is convenient to count modes from this Fermi point and change 
p — > p — pp. Then the vacuum is the state with all non-positive momenta filled by fermions. It is defined as the state 
|0) such that: 

V^|0) = 0, ifp>0, $|0> = 0, if p < 0. (2) 

A state, where in addition to the filled negative modes, m > (— m > 0) consecutive positive (negative) modes are 
also filled (empty) is denoted by |m) (| — m)). 
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Excited states are obtained by a repeated action of fermionic bilinear operators Apgtp^tpq creating particle-hole 
excitations on the vacuum. These operators form an infinite dimensional algebra gl(oo). Coherent states of this 
algebra are the states obtained by the action of a group element on the vacuum 

(m\g = (j n \ e T, p , q A pq :^li> q : _ (3) 

A coherent state represents particle-hole wave-packets. It is characterized by a infinite matrix A pq . The colon in @ 
and throughout the paper denotes the normal ordering with respect to the vacuum : tp^tpq ■ — ip^ipq — (O|-0pV'g|O}- 
Within the algebra of fermionic bilinears one distinguishes two commutative subalgebras. One is generated by 
positive k > (or by negative k < 0) modes of the current 

Jk = ^2 ■ ip^ipp+k [Jk, J-k'] = kdkk'- (4) 
p 

This subalgebra is central in the construction of Ref . 0, Q ■ The other commutative subalgebra will be central in the 
construction of the currentpaper. It is defined in Q13[l. 
An important result |lj, |3( is that the matrix element 

r m (t) = (m\g e~^ k >° tkJ - k \m) (5) 

considered as a function of an infinite number of continuous parameters t = {tk} and of an integer parameter m is the 
tau- function of an integrable hierarchy. r m depends also on the state (g\ (or on the matrix A pq ), but this dependence 
will be omitted in notations throughout the paper. This means that the family of matrix elements r m (t) obeys a 
bilinear identity: 

j e Z k>n (t k -t' k )z\ m (t _ [z]) Tm+n (t / + [z]) z n dz = Q) (6) 

where t ± [z] = {tk ± z~ fc /fc} k = 1,2, . . . and t, t' are two independent sets of parameters. The integral goes along 
a small contour around the origin. 

The bilinear identity © encodes an infinite set of hierarchically structured non-linear differential equations in the 
variables (flows) tk- For a given integer n, this set is called the n-MKP (modified KP) hierarchy. The differential 
equations with respect to the flows tk in a bilinear (Hirota) form are obtained by expanding the integrand of JJjJ in 
powers of z~ l and setting t' k = tk- 

The KP-hierarchy (n = 0) consists of a set of differential equations on r m . The first equation of the KP-hierarchy 
reads 

(Df + 3D 2 2 - 4DiD 3 )t to • r m = 0, (7) 

where Dk is the Hirota derivative in the variable tk defined as a ■ b = (dt k — d t > k ) n a(t)b(t')\ t > k = tk . This equation is 
a bilinear form of the Kadomtsev-Petviashvili equation of plasma physics 0. It gave the name (KP) to the hierarchy. 
The MKP-hierarchy involves a discrete evolution in rn. The first equation of the 1-MKP-hierarchy reads 

(D 2 + D\)T m+l ■ r m = 0. (8) 



B. Evolution in the space of parameters 

Some physical applications arise if tk = ~t*_ k . We write tk — ibtA^k and notice that the operator ^2 k tkJk — 
iSt J A(x)p(x)(dx/2ir) represents the action of the scalar electromagnetic potential A(x) = i^2k^-k elkx 

e -E*t u J h=tr mfA(x)p(x)dx } p( x )=:^(x)xp(x) :, (9) 

which was turned on instantaneously for a short period St. This matrix element is relevant, e.g., for the process when 
an electronic system being initially in an excited state (g\ = (0\g is instantaneously hit by radiation. The tau-function 
then is the amplitude of the probability to find a system in the vacuum (ground) state. 

An alternative interpretation of the formal expression JSJ, common in the literature, emphasizes the commutativity 
of flows with respect to parameters tk- In this interpretation the operator c^- ktkJ ^ k is viewed as an evolution 
operator with a set of "times" tk corresponding to a family of mutually commuting "Hamiltonians" Jk- The excited 
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state (g\ = (m\g evolves in "times". The tau-function is an overlap between the "evolved" excited state (g(t)| and 
the vacuum 

T m =(g(t)\m), {g{t)\ = (m\ge-^ t * J -». (10) 

This terminology, although common, is somewhat confusing. Of course, the "Hamiltonians" , Jk are not Hermitian 
and there is no energy associated with this evolution. The "evolution" occurs in the space of parameters. This may 
be a reason why applications of the powerful mathematical apparatus of soliton theory to electronic physics have been 
limited. 

In the next subsection we introduce the tau-function as a matrix element of the real-time evolution generated by 
the Hermitian Hamiltonian Q. 



C. Evolution in real time 



In physics one is interested in the real-time evolution driven by a Hermitian positive Hamiltonian as Q rather than 
currents 

Can the theory of solitons be used to study the dynamics of fermions in the real space-time, when the evolution is 
described by a unitary operator 

e -iHt +l Px (11) 

corresponding to a positive energy JQ| and momentum 

pel 

In this case the vacuum is selected as a ground state of the Hamiltonian 

H\0) = 0. 

Indeed, the set of Hermitian Hamiltonians 

H h = Y,V k -^l%-, k = 1,2,... (13) 
pel 

generates another commutative subalgebra of fermionic bilinears. Here we denote the set of quantized momenta for 
the system by I, 

Surprisingly, this question has never been addressed systematically. The most relevant papers we know are Refs. 
0,0 where the relation of time dependent Green functions of impenetrable bosons and Green functions of Ising model 
to the Painleve equations has been established. 

In this paper we will show that, indeed, an action of the commutative Hamiltonians 1|13|) generates an integrable 
hierarchy. Specifically, we show that the matrix element 

Tm (t) = (O|0e-£*>o** H * : e -(«+"»M*) . jo), < |a| < 1/2, m is integer (14) 

obeys the MKP hierarchy . 

Namely, the tau- functions r m (t), similarly to JSJ, satisfy a set of equations analogous to |JZJ|SJ|. These equations 
are generated by the identity: 

o T(z + m) 

which is a modified form of ©. The factor T ^(z+m)^ = JX^m^ 1 ^ 2 ^~ ^ replaces z n in © 2 . 



1 Here we take a ^ 0. Although the tau-function obeys the same a-independent equations, it degenerates at a = and can be computed 
by elementary means. 

2 The difference between 1151 and JSJ is inessential insofar as the former may be transformed into the latter by the transformation 
Tm — * ^( m + fc ) j where V(z) = "}2 k t k z k , however, this transformation changes analytical properties of the tau-function as a 
function of the "times" . 
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In the eq. (|14(l the operator ip is defined as a Bose field 

^ = J2i elkxJk+u ° x ( 16 ) 

fe#0 

and 

where ip±(x) = J2±k>o( elkx /k)Jk contain only positive or negative modes respectively. Note that the flow t\ just 
shifts the coordinate x by it\ as e tlHl <p(x)e~ tlHl = tp(x + iti), so we may set x = in some formulas below. In 
physical applications we identify 

time : t — U2, space coordinate in the Galilean frame : x — vpt = —it\, (18) 

where vp = pp/M is a Fermi velocity. 

The operator : e aip ^ x > : in i|14|) is called a boundary condition changing operator. This operator acting on a vacuum 
creates a boundary state 

\B m {x)) =: e -( a + m M*) : |0). (19) 

The meaning of the boundary condition changing operator is clarified by the action of the counting operator n(y) = 
Jq(: ijA (y')ip(y') '■ —(a + m))dy' on the boundary state 

{B m (x)\n(y)\B m (x)) = (a + m) if < x < y, or otherwise. (20) 

In other words, the boundary condition changing operator moves particles towards the point x creating an increase 
of the density p(y) — d y n(y) at that point. 

Boundary states play an important role in 2D-critical phenomena in systems with boundaries. The connection of 
the latter to lD-fermions is following. The fields (e.g. ip(x) and <p(x)) are originally defined on the unit circle. One 
views the circle as a contour embedded in the complex plane z. The fields and their differentials can be analytically 
extended to the exterior of the circle and to its interior. Instead of considering fields defined at the point z' in the 
interior of the circle one views them as defined at the point z = 1/z' of the circle's exterior (Schwarz reflection). From 
this point of view the circle z = e~ lx , x £ R, is seen as a boundary of the holomorphic 2D-field theory defined in its 
exterior. 

In electronic physics boundary states appear in numerous problems related to the Fermi-edge singularities, quantum 
impurity problems and other phenomena related to Orthogonality Catastrophe (|lCj|. see, e.g., We will discuss 
some applications to electronic physics elsewhere |l4j . 

In terms of the boundary state, our tau-function Ijl4(l reads 

T m {t) = {G(t)\B m ), (G(t)\ = {0\ge-^*>° t * H K (21) 

The tau-function depends on two parameters: a fraction a, and an integer m. Below we fix a and consider m as a 
running parameter. 

One can think of this matrix element as follows. An excited coherent state (g\ evolves up to a time t and then is 
measured at a point i by a projection onto a boundary state \B a (x)). 

A comment is in order. The presence of the boundary operator : e aip ^ : changes the momentum quantization. The 
eigenvalues p of the momentum operator ((T2"|> appearing in (|T2l IT^|l are no longer integer, but are integers shifted by 
a: 



peZ + a, (22) 

that is / in 112|) and (|13fl is given by / = Z + a. 

An obvious generalization of this result is a passage from the KP to the 2D Toda hierarchy @. The matrix element 

(B m \e-^k ikHh ge-^* thHk \B m ) 

obeys the 2D- Toda hierarchy with respect to two sets of parameters t}. and i). and a discrete parameter m. Other 
obvious generalization occurs when fermions are placed on a spatial lattice so that the energy spectrum is a periodic 
function of momentum. We do not discuss these generalizations here. 
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The modified bilinear identity i|15|) generates an infinite set of Hirota bilinear equations. It, therefore, obeys KP 
hierarchy with respect to infinitely many flow parameters t^. It also obeys the MKP hierarchy with respect to a 
discrete parameter m characterizing the boundary state. In a special case when only two physical flow parameters 
(|18H t2 and t\ are present we obtain the 1-MKP equation. In the Galilean frame it reads 

{iD t + D 2 X - 2imD x )r m+ i ■ r m = 0. (23) 

The independent variables x and t in this equation are physical space and time coordinates of a ID fcrmionic system. 
Among the hierarchy, this equation seems to be the most important for the study of the dynamics of electronic 
systems. The two functions r m and r m+ i are not independent but are connected by certain analytical properties, thus 
making equation l|23|) closed. In the most interesting situations (corresponding to particular choices of coherent state 
g) the analytical conditions are explicit and correspond to a certain reduction of the MKP-hierarchy. We consider 
coherent states having particle- hole symmetry and the corresponding reduction in Appendix 1X1 

We note here that the matrix element (|14|> give solutions of the MKP hierarchy, which are essentially different from 
the ones given by J5|. For example, we show below that l|14|) has only positive Fourier modes with respect to all 
odd flows iti,it3, . . .. This analytical condition is restrictive. In particular, it does not admit soliton solutions. In 
contrast, it is well known that (J5J can generate soliton solutions (see Sec llllll . 

Some mathematical constructions related of the matrix element (|14|l . although devoted to different problems, have 
appeared recently in papers by A. Orlov j(| and by A. Okounkov, R. Pandharipande, and N. A. Nekrasov Q. 



II. HIROTA EQUATION AND BOSONIZATION 

Before we proceed, we present another form of the bilinear identity i|15|) for a particular case n = 1. Setting 
t = t'°) + [p], t' = t'°) — [q] the factor (z + m)e^<=( tfc ~ tfc ) z! " becomes ^§^j7^h • The integrand has three poles at 
z = 0,p, q,oo and can be computed. It gives the discrete Hirota equation 

(p + m) r m ■ T m+1 ([p] - [q]) - {q + m) T m ([p] - [q]) ■ r m+1 - (p - q) r m ([p]) • r m+ i (-[#]) = 0, (24) 

where we have dropped the term t^ ) in each of the arguments of the r functions. We stress that p and q are shifted 
integers p,geZ + o. 

Hirota's equation, written in the discrete form (|24|l appears as a particular specification of the integral form of the 
bilinear identity ©. In fact, they carry the same information as an integral form of the bilinear identity written for 
an arbitrary n. The latter can be obtained from the former In the rest of this section we will derive the discrete 
form of Hirota's equation (|24|) for tau- functions defined in 114|) . It will be equivalent to proving (|15[) . 

Our main equation 1|23(1 directly follows from Hirota equation written in the form l|24|) . It appears in the leading 
order expansion at large q and p, and setting q = —p. 



A. Bosonization in coordinate space and the tau-function r m =< g(t)\m > 

First, we show how to derive the discrete Hirota equation l|25|) for the the conventional matrix element matrix 
0. Then we derive the discrete Hirota equation l|24|) for our tau-function (|14fl . The first derivation is well known 
P, H, S E3 • We, nevertheless, present it here in order to use it as a guide for the following derivation of Hirota's 
equation for the real time evolution. The comparison of these two derivations is instructive. 

Applying the same trick as is the previous section to © we obtain a discrete Hirota equation for the tau-function 
© 

PT m ■ T m +i([p] - [q]) - qT m ([p] - [q]) ■ r m+ i - (p - q) T m ([p]) • T m+ i (-[?]) = 0. (25) 
To derive it one needs two technical ingredients: bosonization and the Wick's theorem. 

Bosonization formulas 

The Bosonization formulas express the action of fermion operators in terms of bosons. We will need the following 
formulas (briefly derived below) 

^(x)ip(y)\m) = 1— T : e ^~^ : \m), (26) 
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4>(x)\m + 1) = e lx : e'^ x) : \m), <tf(x)\m) = : e"^ : \m + 1). (27) 

Let us now use these bosonization formulas to "fermionize" the tau-function r TO+ i([p] — [q]). We set q — e ix and 
p = e ly , and write the tau-function r m+1 ([p] — [q]) = (m + l\g e~^<= J -*> ; e v(v)-v(x) . | TO _|_ ^ appearing in 
in terms of fermion operators: 



lfljp] - M) = e i(m+1)(:E - a) (l-e 4 ( :E -^)(m + l| 5 (t)V' t (x)V(2/)|m + l}, s(t) = g e - £* (28) 



Wick's theorem 

We apply Wick's theorem to four fermion matrix elements. Let Xi De a linear combination of fermionic modes 
Xi — En r?i an< i IG2) an d (Gi are coherent states (i.e. states of the form ©). The general form of Wick's 
theorem for four fermion matrix elements reads 

(Gi|x t 1 X2xk4|G 2 )(G 1 |G 2 ) = (G 1 | X t 1 X2|G 2 )(G 1 |xk4|G 2 ) + (G 1 | X t 1 X4|G 2 )(G 1 | X2 xi|G 2 ). (29) 

A proof of the Hirota bilinear identity \25\l 

We apply Wick's theorem to the matrix element: 

< m\ip m+1 g(t) ^\x)t/)(y)^} n+1 \rn}(m\g(t)\rn) = (m + %(t)|m + l)(m\g(t) ip^(x)ip(y)\m) 

+ (m + l\g(t)^(x)\m){m\g(t) 4>(y)\m + 1). (30) 

Eq. (|25(l emerges if one writes each of the matrix elements appearing in this identity as r-functions with shifted 
arguments, by further use of the bosonization formulas H26I27[) 3 • 

Derivation of bosonization formulas I26\21\j 

It will be instructive to have a short proof of bosonization formulas (|26I27|) as we are going to use a similar approach 
in the next section. In the following we set for simplicity m = 0. 

Let us project the r.h.s. and the l.h.s. of (|26|l onto a general state with a given number of particles ({x,, y{\\ — 

(0| (j^Li^ {xi)ip{yi)j , consisting of N particle and N hole excitations with coordinates Xi, yi, and compare the 

results. 

First of all, let us use Wick's theorem to calculate the overlap between the state with N particle-hole pairs and the 
vacuum. This overlap is a determinant of a single-particle Green's function (0\^ (x)ip(y)\0) = (1 — e^ x ~ v >)~ x 

.~ N ( 1 \ .^iv n^«,v(e iyi -e^)(e^ - e iXi ) 
({ Xi , Vi } 0) = e^f =I Vi det — = v^Jl<2^ (31) 

where we used the multiplicative form of the Cauchy determinant 

det (^—) = ^<^-^-^\ (32) 
i,j<N \ q t - Pj J \li,j<N W ~ Po ) 

The l.h.s. of (|26|) is computed by extending the previous formula to N + 1 particle-hole pairs and setting xn+i = x 
and y N+ i = y 



(K «}|^(x)^(»)|0> = —— II ( e iy-e^)(e fa -e^) ^ W}|0> ' ^ 



3 This proof appeals to a general view on tau functions as Pliicker coordinates of the Grassmann manifold 2] . 
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The r.h.s. of l|26|) can be computed by (i) noticing that e hv |0) = e^-jO), (ii) pulling the vertex operator e bv ~ to 
the left using the identity 

( p ix piyi \ 
—-—) tfixiWvt) (34) 

and then (iii) using (0|e 6 ^-^ =0. As a result we get 

»}| : e^~^ : |0) = J] L _ e JL _ e J »>l°>" ^ 

i<n y - n -J 

Comparing eqs. I|35|) and l|33|) gives the bosonization formula (|26|) . Eq. (|27|) can be derived in a similar way, or by 
setting x — — ioo in l|26|) and using the fact that ip '(— ioo)\m) = \m + 1). In all formulas x and y carry small positive 
imaginary parts. 



B. Bosonization in momentum space and the tau-function r m =< G(t)|B m > 

Now we show that the matrix element (|14|l obeys the Hirota equation Ij24(l . The steps of the derivation in this 
section follow the steps of the previous section, but technically different. We repeat that in this section all momenta 
are shifted integers p, q E Z + a. 

The important property of the boundary state is displayed by the form of its wave-function in momentum space. 
We have (0|£> m ) = 1 and two important formulas 

(0\^ p \B m ) =/+(p)J_/-( g ), q<a<p, (36) 
p-q 

P±Ht(0\^ p \B m ) = (0\1>ty P \B m+1 ), (37) 
q + m H H 

where 

T(p + to) 



r(m + a)r(p — a) ' 



- r(l — q — to) 

f ™ {q) = T(l-m-a)T(l- q + ay (38) 

The matrix element Ij36|l is also non zero at p = q < a, where it is equal to l 4 . Eq. (|37J) follows from Ij36l38|) with the 
help of the formula 

ti(p)f m (q) ^-^ = f++Mf m+ M- (39) 

Q ~ \~ 777 

We prove these formulas at the end of this section. 

The formula (|36[) is related to a phenomenon referred to as "Orthogonality Catastrophe" 01 ■ Namely, if the length 
of the ring is large, i.e, in the limit p, —q 3> \m\ the formula l|36l) reads 



, , ,+ ,„ i / N,„sin(7ra) 1 f p\ 

n p-q \-qJ 



a+m 

(40) 



The overlap between the excited state and the boundary state vanishes as the momentum of the excited particle 



approaches the Fermi-point (O|')/'J , 0p|J5 m ) ~ p a+rn as p — » 0. 



4 The formula 1361 is understood as being analytically continued from the real axis to the complex plane of p and q. In particular, at 
p = q it should be understood as a limit p — * q at fixed q. 
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Consider now a generic excited state ({pi, qi}\ = (0 Yii=i ^q^Pi characterized by a set of iV-particles with momenta 
Pi > a and N holes with momenta qi < a. The overlap of this excited state with the boundary state is obtained from 
the Wick's theorem and ea. (|36|l 

{{piAi}\B m ) = det , m\ iP Pi \B m ) = .det ( f+( Pi )—!—f-( qj 

i,3<N H3 i,J<N \ pi — Qj 

= ii fMfM "'tV'iv^ < 4i » 

where in the last line we again used the Cauchy determinant formula l|32|) . 

Let us note that up to the multiplicative factor Yli fm(Pi)fm(li)^ the overlap (|41|l between the excited state and the 
boundary state written in momentum space is given by the Cauchy determinant, i.e., has the same form as the similar 
overlap written in real space The determinants become identical under a substitution pi — e lXi and qi = e lVi . 

We observe an interesting phenomenon: the boundary state \B m > interchanges momentum and coordinate spaces. 
The transmutation 

coordinate space <-> momentum space, (42) 
Fermi sea <-> Boundary state (43) 

is the reason why different matrix elements (JSJ I14|) obey similar non-linear equations. There is no one-to-one corre- 
spondence between coordinate and momentum space - the space is continuous and finite while the momentum space 
is discrete and infinite. 

Following the logic of the previous section we proceed in two steps. First, we establish the bosonization formulas in 
momentum space. They will allow us to rewrite tau-functions in eq. (|24|l as matrix elements of fermionic operators. 
Then eq. I|24|) is equivalent to a particular form of the Wick's theorem. 

Bosonization formulas 

Let us introduce a bosonic field associated with the set of commuting Hamiltonians l|13fl 

= E S -y H * = E ln d - - s ) ■ 4>ll> P ■ ■ (44) 

k>0 P 

Bosonization in the momentum space is summarized by formulas 

4*l> P \B m ) = ^^ q \ ^)-^)\ Bm ), (45) 

and 

km (p + m)V'JV'p|-B m ) = (-ir^^/- +1 («/)e*^| J B Tn+1 >. (46) 

p — m 7T 

These formulas will be derived below. We remark here that (|45I46|I should be understood as analytically continued 
formulas for matrix elements. For example, one should first calculate the overlap of the l.h.s. of l|4()|) with some 
bra-state, then analytically continue it from p G Z + a to the complex p-plane and only then take a limit p — > — m. 
Once bosonization formulas in momentum space (|45I46|I are established we are ready to prove the Hirota equation 

pi|i. 

Wick's theorem 

First, we apply the Wick theorem to the 4-fermion matrix element 
(G(t)\^ p ^ P \B m )(G(t)\B m ) = {G{t)\^ p \B m ){G{t)\^y P \B m ) - (G(t)\^ P \B m ){G(t)\^ p \B m ). (47) 



Next we multiply both sides by p +yv } send P — > — m, and use eq. (]46fl. We get 
^(G(t)| e *W)^^|i? m+ i}(G(t)|B m ) 
= <G(t)|^^|S m )<G(t)|e*^)|S m+1 ) - ^±iM (G(t )|e*(«)|S m+1 )(G(t)|^V P |S m ). (48) 

/m+iw) 

The next step is to use the bosonization formula (|45|l 

«^£±^ (G( t)| e *W) e *(9)-*(p)|B Tn+1 )( G (t)|S Tn ) 
p - Q q + m 

= (G(t)|e*W-*W|S m )(G(t)|e*W)| J B m+1 ) - E^± {G (t)\e^\B m+1 )(G(t)\e^-^\B m ). (49) 

g + m p — Q 

We recognize the matrix elements in ()49|l as tau-functions l|14ll with shifted arguments, e.g., (G(t) e*^ \B m+ i) — 
r m+ i(— [q]). The last step is to take a limit Q — > oo. In this limit $(Q) — > and 149|) brings us to (|24|l . which 
completes the proof. 

Derivation of formulas for matrix elements lcSbl38{) 
Now we prove the formulae i|3t)l38|) . First, we find using (|34fl that in coordinate representation 

/ -| ix \ m ~\~ a iy / -| ix \ 

(0|^(^(y)|i3 m ) = (0|^(^My)|0) — ^- =— ^— 7^7 • (5°) 



l_ e iy e iy _ e ix li _ g ly 



We consider (p — q)(0|'(/'J'(/ ; pl^m) and pass to the coordinate representation (p — q){0\ipj l ip p \B n 
(id x + id y )(0\^ ' (x)ijj{y)\B m ) . Then we calculate 



(id x + id v ) 



^iy / gj.-r \ rri-\-a 



e W — e ix V 1 — g}V 



= {m + a)(l - e* B ) m+a - 1 e iv (l - e i v)-( a + m + 1 ) . 



We observe that the right hand size factorizes. Taking an inverse Fourier transform we obtain (|36[) with 

/-(?)= /' 27r e ^(i_ e « r +a-i^ ) y+(p) = (m + o) / ,2 ' T e -(p-i)v(i_ e ^)-(™+"+i)^. 
Jo 27r Jo 2?r 

The calculation of these integrals gives for . These formulas can also be found in . 

Derivation of bosonization formulas \4^4®il 
To prove the formula 1|45(1 we evaluate the projection of its r.h.s. on the state ({pi, qi}\. With the help of the formula 

e»W^ e -»W = ^£z£y^ (51) 
analogous to l|34|) and using (0|$(s) = we obtain 

({ Pi , qi }\e^-^\B m ) = ({ Pi , qi }\B m ) J] (P-^fa-g') (52) 

On the other hand, using the multiplicative representation for the l.h.s. of l|45|) similar to l|41l) we obtain 

({ Pi , qi }\^ p \B m ) = . det C/+( Pi )^_/-( ft ) > ) = ({ Pt ,q t }\B m ) mp)fm{q) J] (P-P'ffi-gi) (53) 
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where pn+i = p and qN+i = q- Comparing (|53|) and i|52|) we prove l|45|) . 
To obtain (|46[) . we first use i|39|) to obtain 



({ Pi , qi }\B m+1 ) = ({p<,q t }\B m ) I] ?*±^. (54) 

i<N q * + m 

Then we use (|54|) and l|53() to proceed as follows 

= ( {ft>ft} |^|^ > ^__JI— — . (55) 

The r.h.s. of (|55|1 effectively does not depend on p and we take its limit as p — » — m. The product in the last line 
disappears and we conclude that 

e^\B m+1 ) = lim -jL-^jB^). (56 ) 

Jm(Q) m Jrn(P) 



Finally, using 

lim (p + TO )/+(p) = (-ir+ 1 (m + a)^^, ,r,7) 

p — 171 

we arrive to l|4fijl. 



III. MULTI-PERIODIC SOLUTIONS 



In this section we present explicit formulas for multi-phase solutions evolving in real time. These formulas reveal 
the structure of the tau-function and provide another direct proof of Hirota's equation. 

A set of excited states \{pi, q{\) characterized by momenta of N particles pi > a and holes qi < a, i = 1, . . . N form 
an orthonormal basis for the Fock space of fermions. In this basis the Hamiltonians (|13|1 are diagonal: Hk\{pi, qi}) — 
J2i(Pi ~ Qi)\{Pii Ii})- Using this basis, we decompose 

T ">=E E e~ Ek ' itkiP *~ qt H9\{Pi,cii}){{Pi,Qi}\B m ). (58) 

N>0 {p,,?,} 

The matrix elements entering this expression are determinants of two-particle matrix elements. The first one is 

{g\{pi,<n}) = 4 e Mpi,gj 5 A p, q = (g\p,q)- (59) 

The second one is the Schur function given by eq. (|41|l 

({Pi, qi }\B n ) = det (—!—) J] fMfM- (60) 
i,j<N \pi - q } ' i± 



Introducing 



where 



we write the tau-function in the form 



( ({p i , g 0, 5 ) = det(4™).) ) 

42 = fnMA P , q f m (q) 



V V e-EM'*-*),,,^,,},^ _!_ . (61) 
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It is obvious from this form and inequality Pi > a > qi that the tau function has only positive Fourier components 
with respect to all odd flows iti,it 2 , ■ ■ ■■ Indeed, for all odd k we haver p\ — qf > 0. The bilinear identity (|T3|l follows 
directly from this representation if one uses identities for Schur functions and the property 

We do not present this proof of l|24|l here 5 . 

If the matrix A pq defining the state (g\ by l|59|l is a finite matrix, the sum in H61|) truncates. It represents a 
multi-phase solution. In this case the determinant formula holds 

r TO (t,o) = det (Sa + f+(pi)e- e ^^^e^f-( qj )) , e(p) = ^f^p*. (63) 
Here Ai = A Piqi . Note that a similar determinant representation holds for JSJ). Namely, 

r m (t) = det (Sa +p^ e -'to)-^—e<Mq? +1 ) , e(p) = £t*p*. (64) 

i,3<N \ Pi — Qj / ^ 

Although l|63|l and (|64|l are very similar there is an essential difference between them. Namely, in (|63|l the parameters 
Pi and qi are momenta of particles and holes and are subject to "Fermi sea" restriction Pi > a > qi . On the other 
hand, in (|64|l p = e lx and q = e ty represent complex coordinates of particles and holes and are not subject to this 
restriction. Due to the restriction on the values of pi and qi the tau function (|63|) has only positive Fourier components 
in all odd flows. 

The simplest 1-phase and 2-phase solutions are written below as a function of space t\ = ix and time t% = —it. 

- 1-phase solution describes one particle-hole excitation: 

r m = 1 + A ^ P ^ m ^ e l{p - q]{x ' vt) . (65) 
p-q 

The solution is characterized by two positive numbers p and —q from Z + a labelling particle-hole momenta and 
an arbitrary m-independent constant A, which may depend on p, q. It propagates with velocity v = p + q. 

- 2-phase solution (two particle-hole excitations): 

T m = 1 I yl /"( Pl )^( gl ) c '(Pi-<n)*-»(p?-q?)t | A Jm(P2)fm(q2) c i( P2 - q2 ) x -i (l Z- q Z)t /ggN 

Pi - qi l>i 'i> 

+ AiA 2 det ( fmil^fmiPi) \ c i(pi-Q 1 )(x-v 1 t) c i(p 2 -Q2)(x-v 2 t) _ 

i,j=l,2 \ pi - qj J 

The multi-phase solutions in (|64|l allow for rational limit. Keeping Vi = pi + qi constant while taking the limit 
Pi — qi — > 0, one obtains rational solutions corresponding to iV-solitons with u, being asymptotic velocities of solitons. 
However, in (|63|l and its particular cases I|65I66|I this limit does not exist. Indeed, under the condition Pi > a > qi 
the limit pi — qi — > implies Pi,qt — > and therefore it is impossible to have finite Vi in this limit. We see that the 
matrix element <|14[) does not contain soliton solutions of the MKP hierarchy. 

IV. CONCLUSION 

We have shown that the real time dynamics of certain matrix elements (|14f) of one-dimensional (non-relativistic) 
fermions obeys non- linear integrable equations of the modified KP hierarchy. This evolution is driven by a commutative 
set of Hermitian Hamiltonians Ijl3(l 



H, = 



( d \ k dx 



5 An observation that the series 1611 for 1141 obeys KP-hierarchy, i.e., eq. 1151 at n = has been made in Ref.|(j. 
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It has to be contrasted with a traditional representation of the tau-function of the modified KP - hierarchy (JSJ) 
P, Q in terms of fermions. The latter describes an evolution driven by a commutative set of (non-Hermitian) current 
operators |@J 

Jk = j e lkx : i>Hx)iP(x) : 

In these two cases the flows have different meaning. The former are coordinates of physical space-time, while the 
latter are deformation parameters of the coherent state. 

The matrix element we studied in this paper is understood as an overlap of the evolved coherent fermionic state 
(G(t)| with a "probe" boundary state \B m ). This matrix element appears in variety of physical problems involving 
the phenomena referred to as the Orthogonality Catastrophe. Among them are problems related with the Fermi-edge 
singularity, quantum impurity and tunneling in effectively one dimensional systems. 

The fact that questions of quantum dynamics in one dimension are linked to the tau-function, indicates that the 
kinetics are essentially non-linear and are a subject to instabilities inherited from the non-linear dynamics. In a 
forthcoming paper we describe the Whitham approach to shock wave type solutions for electronic systems. 

An essential element of our proof of the MKP hierarchy is the representation of a fermionic mode by a Bose field 
(I44I46|I in the action on a boundary state, 

This is a representation in momentum space. It is valid in a subspace of the Fock space consisting of excitations on 
top of the boundary state. 

Again it has to be contrasted with the traditional "bosonization" - representation of Fermi-field in terms of a 
Bose-field with respect to the Fermi- vacuum 

This is a representation in coordinate space. It is valid in a subspace of the Fock space consisting of excitations on 
top of the Fermi- vacuum. 

We emphasize that although the matrix elements iJSJ and (|14|) obey similar MKP equations they give essentially 
different solutions of those equations. For example, (|14|) gives only solutions satisfying analyticity requirements in 
odd flows. The latter condition excludes soliton solutions. 

Finally, we emphasize that the Schur decomposition H61(l of the tau-functions similar to (|14jl (an essential element 
of the proof) has appeared in in the problem of counting of Hurwitz numbers and in 01 in constructions of tau- 
functions of hypergeometric type. Similar objects arise in supersymmetric gauge theories 0, 1 1 5| and also in matrix 
quantum mechanics. 
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APPENDIX A: REDUCTION FROM THE PARTICLE-HOLE SYMMETRY 

Let us treat the tau-function as a function of space-time l|18|) . At first glance, the equation lll'-'il) is not closed. It 
is written on two functions r m and r m+ i. In fact, both functions are determined by the initial state (g\ and are not 
independent. The relation between these functions determined by formulas (|61I62|I makes l|23|) a closed equation. 

However, in the case of generic g the relations between r m and r m+ i are difficult to formulate. The situation is 
simplified when the initial state {g\ has some special form, or is restricted by symmetries. One may think about these 
cases as reductions of MKP hierarchy. 

The multi-phase solution l|tjl[) is the simplest (finite-dimensional) reduction. In this case the tau-function, treated 
as a function of the variable z = e~ lx is a polynomial in z . 
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A more interesting reduction having direct physical applications occurs when the parameter a + m takes a "sym- 
metric" value such that a + m = — (a + m + 1). We set 6 

a = -1/2, m = 0. 

In addition we require the state (g\ to be invariant under a charge conjugation 

e^0( P ^) ^ |p, q) = -{g\- q , -p), (Al) 

where xq is some real constant. This condition reflects a a particle-hole symmetry: ?/>J — > "0-?: — * ^-y 

Under these conditions the tau-function treated as a function of real space-time coordinates (|18|l has a reflection 
property 



t (x) =n(-x + 2x ), (A2) 

where r stands for complex conjugation of the function. This property makes the equation (|23|l closed with the second 
complex function T\ being completely determined by t . 

The analytical meaning of the condition l|A2|) becomes clear if one considers an extension of the tau-function to the 
complex plane z from the unit circle z — e~ l ( x ~ x °K Then the reflection condition mea ns that the coefficients of the 
formal series l|61[l for To and t\ are complex conjugated (we use a notation f(z) = f(z)) 

T (z)=n(z), (A3) 

This follows from the property 

/i + (p)/r(?) = / + (-?)/o"(-p) 

held at a = -1/2. 

The holomorphic tau-function t(z) = t (z) is defined in the exterior of the unit disk, while its Schwarz reflection 
t(z) = f 1(1/2:) s defined in the interior of the disk. The reflection condition l|A3(l means that the holomorphic 
tau-function, and its Schwarz reflection, can be glued along the unit circle forming a function analytical in some 
neighborhood of the circle. 

Particle-hole symmetry IjAljl occurs when A pq is a Toeplitz matrix, i.e., it depends on the difference p — q. In this 
case the state (g\ is given by 

g = ^jA(x) P (x)dx^ A(x-xo) = -A(-x + xo), (A4) 
where A(x) - the symbol of the Toeplitz matrix i.e., a Fourier transform of A pq . It is an odd real function. 
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